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ABSTKACT 

Let  a  Iteas  be  put  on  test  st  the  outset,  and  suppose  an  item  Is  not 
replaced  upon  failure.  Assume  an  exponential  failure  dlatrlbutlon 
Fg(t)  -  1  -  exp(-t/6).  A  tlae  truncated  aequentlal  procedure  for  testing 
B ^ versus  6  s 0^  le  developed.  This  procedure  allows  a  quick 
rejection  of  when  Is  true,  but  provides  an  accurate  Interval  eatlaate 
of  6  when  Is  accepted  after  the  test  has  been  established. 

AMS  Subject  Classification;  Prlaary  62L10 

Key  Words  and  Phrases;  Tlae  truncation;  Gauaslan  process;  Wiener  measure; 

Weak  convergence;  Stopping  rule;  Confidence  Interval 
Average  aaapllng  tlae. 
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1.  INTRODUCTKW 


In  the  quality  control  of  industrial  production,  it  Is  a  coonon  practice 
to  test  hypotheses  to  make  a  decision  about  acceptance  or  rejection  of  a  pro¬ 
duction  batch.  Generally,  a  sequential  testing  procedure  requires  less  data 
to  reach  a  decision  than  a  fixed  sample  or  fixed  time  testing  procedure.  Also, 
in  many  situations,  the  experimenter  may  want  to  obtain  a  point  or  an  interval 
estimate  of  some  parameter  after  a  test  of  hypothesis  has  been  established. 

(for  example,  see  Slegmund  (1977,1978)). 

In  this  paper,  ve  will  consider  a  commonly  used  life  testing  procedure 
In  Industrial  settings.  Let  m  Items  be  Independently  put  on  test  at  the  out¬ 
set,  and  when  an  Item  falls.  It  Is  not  replaced.  Suppose  that  the  failure 
time  Is  exponentially  distributed  trlth  distribution  function  F.(t) ~ 
l-exp(-t/e),  t 2  0.  The  hypotheses  which  are  to  be  tested  are 

(1.1)  Hq:  0^00  versus  Hj^;  6  se^,  where  ej^<  8^. 

We  assume  that  for  certain  reasons  this  life  test  must  be  terminated  at 
time  t^  If  not  all  of  the  Items  have  failed  by  tQ.  (This  type  of  censor¬ 
ship  Is  often  called  Type  I  censoring.)  The  time  t^  can  be  preassigned 
or  determined  by  the  specified  probabilities  a  and  8  for  testing  (1.1). 

The  maximum  likelihood  estimator  of  6  under  the  above  type  of  censor¬ 
ship  with  censoring  time  t(^tg)  was  given  by  Bartlett  (1953)  as 

m 

I  V,  ■¥  (m-n  {t))t 

6  (t)  -  ^ - , 

N  (t) 

B 

where  N  (t)  denotes  the  number  of  failures  at  or  before  time  t  and  the 

B 


2 


U^'s  denote  the  observed  failure  tlaes  In  Increasing  order.  For  convenience, 
let  us  define  6^(c)  ■  **  if  N^(t)  •  0.  Bartholoaew  (1963)  obtained  the  con- 

4k 

ditional  probability  that  6^(t)  2  c  given  N^(t)  >  0,  when  c  >  0.  Yang  and 
Sirvanci  (1977)  showed  that  the  eatinwtor  0^(t)  is  consistent  and,  %rtien 
properly  standardised,  is  asyaptotically  noraal  for  each  fixed  t  s  t^. 
Recently,  Spurrier  and  Wei  (1979)  have  used  9^(1^)  as  the  test  statistic 
for  testing  (1.1)  and  denonstrated  certain  advantages  over  Epstein's  procedure 
(1934)  and  the  pure  sequential  probability  ratio  test  given  in  Ghosh  (1970, 
pp.  193-96).  Note  that  Spurrier  and  Wei's  test  is  a  test  based  on  the  fixed 
length  of  tiac  t^. 

In  this  article,  a  one-sided  sequential  test  procedure  based  on  the 
random  function  6^  is  presented  for  testing  (1.1).  This  procedure  allows 

e 

a  quick  rejection  of  %ihen  is  true,  but  provides  an  accurate  interval 
estimator  of  6  when  is  accepted.  Siegmund  (1977)  studied  the  same  problem 
for  different  expcrisMntal  settings.  In  Section  2,  we  prove  that  the 
standardized  process  6  converges  in  distribution  to  a  Gaussian  process  as 
m  -*'  <».  The  large  sample  approximation  of  the  power  function  and  the  average 
sampling  time  (AST)  are  presented  in  Sections  3  and  4,  respectively.  An 
Interval  estimate  of  6  after  testing  is  given  in  Section  S.  A  simple 
example  is  provided  for  illustration  purposes  in  Section  6. 

2.  LARGE  SAMPLE  APPROXINATIW  OF  6  BY  A 
GAUSSIAN  PROCESS 

To  test  the  hypothesis  (1.1),  we  define  the  following  stopping  rule: 

A 

Let  T  be  the  smallest  value  of  t  such  that  9^(f)  ^  Given  t^,  stop 

testing  at  nin(tQ,T)  and  reject  if  and  only  if  T  s  t^. 


We  now  prove  that  the  distribution  of  the  standardised  9  converges 


weakly  to  a  Gaussian  measure 


Theorem  1.  Consider  6  (t)  to  be  a  function  of  t  c  (0,t>].  Then 


m'(6  -  6)  converges  in  distribution  to  a  Gaussian  process  Y  as  m 


where  EY(t)  -  0  and  cov(Y(s),Y(t))  •  6  /F.(max(s,t)),  s.t  c  (O.t^] 


Proof.  The  total  time  on  test  at  time  t  Is  given  by 

N^(t) 

♦.(t)  -  I  U.  ♦  (m-N,(t))t 

■  1.1  ^  ■ 


First,  we  prove  that  the  random  element  C 


verges  in  distribution  to  a  Gaussian  process  Z  with  EZ(t)  ■  0  and 
2 

cov(Z(s)  ,Z(t))  ■  6  F.  (mln(s,t))  as  m  -»  •>,  where  s,t  c  (See 


Billingsley  (19h8)  for  notation.) 


^  s.  &  t-,  where  k  Is  an  arbitrary  positive  Integer 


X  be  the  hypothetlcel  data  from  the  exponential 


distribution  with  mean  9.  Note  that  because  of  the  time  truncation  t.  not 


every  X  can  be  observed.  Define 


(X,-(8.+«))I 


%fhere  I,  denotes  the  Indicator  function,  1  ■  1 


It  can  be  shown  that  EV 


-  e^(Fg(s  )-Fg(Sj_j)),  J  d  1,  and  E(wJ*'wJ^  ')-0,  t  d  t' .  By  the  Central 

Limit  Theorem,  the  k-dlmenslonal  random  vector  m  L  )i*  converges  In 

1-1 

distribution  to  a  multivariate  normal  random  vector  with  mean  £  and 

,  where  W.  -  . end  t 


covariance 


dlag(o,  ,0 


denoting  matrix  transpose 
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j  ^ 

ut  -  I  I  -  ♦«(».)  -  ew  (sj  «nd  V  -  (V 
■  Z-l  1-1  ^  “  J  ■  J  ^  ' 


(1)  v(k)v. 

9  •••  9  *_  /  • 

n  iB 


Then  m  V  converges  in  distribution  to  s  nultlvsrlste  nomsl  rsndon  vector 

iS 

with  oMcn  0  and  covariance  aatrlx  |  *  (*44)1,  v  **•  vhere 

1  Ij  ll*K 


•ij  “ 


e  f^(8j)  .  j  s  1 
®  ^0^*1^  *  J  ^  ^'■•3  ■  If.tk. 


Let  0  5  <  r  <  s  t^.  Then  C^(r)  -  -  m"’  J  and 

1-1 

C,(rj)  -  t,(r)  -  .her.  -  (X,  -  (r*«)>It,^  +  (r-tj) 

-<'-'l>'cx^<rj]  *1  ■  “’’'[riXj.rj]  *  <V^>  '  <V'>*tX,  <  r)' 

It  la  straightforward  to  show  that  Et^  -  -  0,  -  0, 

2  2 

var  ■  ®  (e*p(-rj^/0)-e*p(-r/e)X  var^j^  -  6  (exp(-r/e)-exp(-r2/8)) ,  and 
Et^8^  -  (r-r2)^e^(exp(-r/e)-«xp(-r2/6)).  Mow, 


Ef(C  (r)-€.(r 

m  m 


,))‘‘(C  (OKCr))^]  -  m"^E[(  I  xYi  I  iYl 

1-1  J-1  J 


-l_r  2.2,  .  -1,  -,_r  2.2, 

-  B  ECtj^j]  +  a  (b-1)E[t2823 


-  B“^(r-rj^)^e^(exp(-r/0)-exp(-r2/8))  +  a“^(B-l)var  •  var82 
s  a  ^(r2-r2)^8^(exp(-r2/8)-exp(-r2/8))  +  B~^(B-l)(exp(-rj/e) 

-exp(-r2/8))^e^ 

s  8*[(r2-r2)^/8^  +  (exp(-r2/8)-e*p(-r2/8))^] 
s  e*[(r2/8-rj/8)  +  (exp(-r2/8)-exp(-r2/8))]^ 

-  8*[(Fg(r2)+r2/8)  -  (Fg(r2)+r2/0) ]^. 


Thus,  by  Theorem  15.6  of  Billingsley  (1968),  the  randoa  eleaent  ( 

M 

convergea  In  distribution  to  a  Gaussian  process  Z  with  EZ(t)  -  0  and 

9  .1 

cov(Z(s),Z(t))  8  F  (ain(s,t))  as  a  Since  the  rendoa  element  mN  converges 

6 

In  probability  to  a  constant  function  l/F^  In  D(0,tQ3,  by  Theorem  4.4  and 


Corollary  1  of  Theorea  S.l  of  Bllllngaley  (1968),  the  raodon  elenent 


C  *  -  6)  converges  In  distribution  to  a  Gauaalan  process  Y 


with  EY(t)  ■  0  and  cov(Y(s) ,Y(t))  ■  6  /F. (aax(s,t))  as  a 


Theorea  1  will  be  used  to  obtain  the  approximate  power  function 


for  the  test  of  (1.1)  for  large  m  in  che  next  section,  to  calculate 


the  average  sampling  tine  in  Section  4,  and  to  obtain  Interval  estiaates 


of  0  after  testing  in  Section  5.  First,  a  standard  kind  of  tlae 


transforaation  for  the  process  Z  in  the  proof  of  Theorem  1  is  needed 


1  for  future  reference 


1.  Let  Z  be  a  Gaussian  process  with  EZ(t)  >  C  and  cov(Z(s) ,Z(t)) 


s  (.  (0,F.  (t^)].  Then  p(s)  is  a  standardised 


is  the  Inverse  function  of  F, 


Wiener  process 


3.  ASYMPTOTIC  APPROXIMATION  OF  POWER  FUNCTION 


Since  N  (0)  ■  0,  by  our  convention  in  Section  1,  0  (0)  is  defined 


0.  Therefore,  we  define  Z(0)  ■  b,  where  b  is  an 


arbitrary  positive  number  which  can  be  determined  from  the  specified  error 
probabilities  a  and  $  mb  a  result  of  Theorea  2. 


P.  (accept  H_)  -  Pft(T>  t„) 


6 


u 

where  c*(0)  ■  ■  (c  -  0)  and  ♦  is  the  distribution  function  of  the  standard 

B  S 

nomal  distribution. 


Proof.  By  Theorem  1,  we  have  for  the  Gaussian  process  Z  of  Lwaa  1, 
Pg (accept  Hq)  -  -  0)  >  c^(e),  for  all  t  c  CO.tQ]) 

m>»  P.  (1-^7^  >  c'(0),  for  all  t  t  [0,t-3)  trhlch  by  Leans  1  becomes 
0  Fgit/  m  u 

p  (Misi  >  c’(0)  ,  for  all  s  c  [0,F-(t-)]) 
o  8  m  o  u 

■  Pg(p(8)  >  c^(O)8/0,  for  all  s  c  [Q.F^Ctp)]) 
b-c* (0)8 

-  Pg(p(8)  <  - 1 -  ,  for  all  8  €  [0,Fg(tQ)]) 


where  the  last  equality  is  obtained  from  an  application  of  the  results  on 
page  348  of  Shepp  (1966).  /// 

Therefore,  in  order  to  determine  the  constants  c  and  b  required  for 

ID 

testing  (1.1)  with  the  specified  error  probabilities  a  and  8,  we  must  solve 


the  following  two  nonlinear  equations: 


-  “^>(2  ‘:;(6o)b/eo^) 


‘  * \ 


A  solution  to  (3.1)  tuy  be  easily  found  by  numerical  analysis  techniques. 
An  example  is  given  in  Section  6. 
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4.  ASYMPTOTIC  APPROXIMATION  OF  DISTRIBUTION  OF  T 


The  approximate  distribution  of  the  tine  T  for  large  m  may  be  readily 
obtained  from  the  results  of  Sections  2  and  3.  The  distribution  function 
of  T  is  given  in  Theorem  3. 


Theorem  3.  When  m  is  large. 


C,(t)  . 


where  t  >  0. 


/cI(fi)Ffl(t)  *b/e\ 


>c^(e)Fe(t)-b/e' 


Proof.  From  Theorem  1,  T.«nnia  2,  and  equation  (17.1)  of  Shepp  (1966), 


Pg(Tit)  -  1  -  Pg(T>t)  -  1  -  Pg  (*^(e^(u)-0)  >  c^(e)  for  all  u€[0,t]) 

-  1  -  >  c;(e),  for  all  uc  C0,t]) 

6 

.  1  _  p  (  Qpix).  >  c’(e),  for  all  V  c  C0,F^(t)]) 

6  V  m  0 


1  -  Pg(p(v)  < 


b-c'  (e)v 


,  for  all  ve  [0,Fg(t)]) 


1  -  {♦ 


^b-c;(e)F,(t)^  /2c;(e)b\  /■  c;(e)P,aHb\ 

(  e<P.(t))'*/  9^  /  (  9(F.(t))'’  ) 


Note  that  T  may  not  be  a  proper  random  variable,  i.e.  Fg(T  •  •)  can 
be  positive.  When  c'(6)>0,  G-(t)  is  very  similar  to  the  inverse 

n  V 

Gaussian  distribution  function  (Shuster  (1968)). 

The  average  sampling  tine  (AST)  for  our  sequential  testing  procedure 
is  given  by 
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f  0 

(4.1)  I  tdGg(t)  -Upa  -  GgCtj,)). 

The  integral  In  (4.1)  cannot  be  obtained  in  closed  form,  but  is  easily 
integrated  computationally  by  the  Gaussian  or  other  quadratures. 

5.  THE  INTERVAL  ESTIMATION  OF  0  AFTER  TESTING 

For  the  completeness  of  this  paper,  we  utilize  Siegmund's  (1978)  method 
of  obtaining  a  confidence  interval  of  a  parameter  8  for  an  arbitrary  parent 
distribution  after  a  test  of  hypothesis  of  6  has  been  established.  First, 
points  on  the  stopping  boundary  are  ordered  in  a  counterclockwise  direction. 

Then  a  lower  (1  -  y)  confidence  bound  for  0  is  the  smallest  value  of  6 
which  gives  probability  at  least  y  to  the  event  that  the  test  terminates  at 
a  boundary  point  at  least  as  large  as  that  actually  observed  in  the  ordering. 

An  upper  confidence  bound  is  defined  similarly. 

For  our  situation,  an  interval  estimate  of  6  when  is  accepted  is  much 
more  desirable  in  practice  than  when  is  accepted.  When  is  accepted,  further 
development  and  testing  irtll  be  necessary  so  that  it  is  more  Important  to 
reach  a  rejection  decision  as  soon  as  possible  than  to  give  an  accurate 
estimate  of  0.  The  following  theorem  gives  a  (1  -  2y)  confidence  interval 
for  0  after  Hq  Is  accepted.  The  proof  of  this  theorem  is  straightforward. 

A 

Theorem  4.  When  Hq  Is  accepted,  let  the  observed  value  of 
Then  a  lower  (1  -  y)  confidence  bound  is  given  by 

A 

(5.1)  0(d)  -  inf  {0:  Pg(T  2  t^  and  2  d)  i  y)  ; 


and  an  upper  (1  >  y)  confidence  bound  is  given  by 
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(5.2)  e(d)  -  8up{e:  [Pg(T<tQ)  +  (T  2  and  6^(1^)  <  d)  ]  2  y)  • 

Finding  the  upper  and  lower  bounds  for  6  involves  computing  the  three 

probabilities  in  (5.1)  and  (5.2).  Since  Po(T2t-  and  0  (t-)  2  d) 

0  0  m  0 

A  A 

“  ^  ^  ^0^^^  ^0  ^d),the  only  quantity 

which  we  are  unable  to  compute  for  large  m  from  results  in  the  previous 

(T  *  tg  and  ®jjj(^q)  ^  d) .  The  following  corollary  to  previous 
results  gives  an  approximate  expression  for  this  probability. 


Corollary  1.  When  m  Is  large. 


Pg(T<tg  and  e^(tQ)2d). 


where  d^  (0)-m**(d-e  ).g(t)  - 
in  Section  2. 


(e)r^<t)-d;(e)F^(t„)tb\ 


g(t)dt. 


,  G  is  defined  In  Theorem  3,  and  b  Is  defined 


Proof.  Write  Pe(0„(tg)2^d  and  T  <  tg)  -  Pg  ^®m  V  ^ ‘*1  ^  ^  ^  V 

•^2  ^^T<tg) 

tg 

Now,  letting  s  “  Fg (t) ,  and  Sg  -  Fg (tg) ,  we  obtain 

ep(s-) 

1  "sg  I  ®P<®^  "  ®‘=m^®^^ 

-  Pg(p(8g)  2  Sgd^  (e)/e  I  p(s)  -  c^(e)s/e) 

"  ^  (so‘^,Jl(®>'Si<®)s)/e  I  p(0)-b/e) 
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■*( 


b-s  d' (0)+c'  (e)s 
On  in 


O/Sq-s 


where  we  have  used  Che  stationary  increments  property  of  p(s).  /// 

The  Integral  in  Corollary  1  can  be  easily  computed  by  numerical 
integration  procedures  in  practice,  thus  giving  approximate  values  for 
(5.1)  and  (5.2)  for  large  m. 


6.  AN  EXAMPLE 

In  this  section,  we  give  a  simple  example  of  computing  the  initial 

value  c  ,b,  and  the  AST  for  illustration  purposes.  We  let  m  -  100, 
m 

6^  -  1.5,  “  1.0,  a  ■  .05,  and  g  ■  .1.  To  attain  given  values 

a  -  P*  (reject  H^,)  and  B  “  (reject  H.),  the  pair  of  nonlinear  equations 

Go  u  1 

in  (3.1)  are  solved.  If  t^  is  preassigned  t^  »  2.079  such  that  Fe^(tQ)  *  .75, 

then  c  ■  1.2A69  and  b  •  1.2436.  We  also  report  the  test  based  on  0  (t*) 
m  m 

for  a  fixed  length  of  time  t*  (Spurrier  and  Wei  (1979))  for  comparison 
purposes.  For  the  same  a  and  g  values,  t*  *  1.0942.  Table  1  gives  the 
average  sampling  times  (AST)  of  our  sequential  testing  procedure  for  several 
G  values,  computed  by  Che  24-polnC  Gaussian  quadrature  formulas.  Under 
the  alternative  hypothesis,  the  AST’s  are  considerably  smaller  than  t* 
for  the  fixed  tine  test.  Under  the  AST  is  larger  than  t*.  However, 
this  is  actually  an  advantage  because  we  would  like  to  have  an  accurate 
estimate  of  G  when  is  accepted. 
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Table  1.  THE  AVERAGE  SAMPLING  TIME 


Iq  -  2.079  . 

t*  > 

1.094 

e 

0.2 

0.4 

0.6 

0.8 

1.0 

1.5 

AST 

0.017 

0.060 

0.129 

0.269 

0.779 

2.009 

7.  REMARKS 

In  practice,  we  perform  the  sequential  procedure  proposed  in  this 

A 

article  by  computing  0  (t)  periodically  with  small  Increments  of  time.  By 

m 

doing  this  our  test  can  be  treated  as  an  approximation  to  repeated  significance 
teats  with  fixed  length  of  time. 

When  t^  is  not  predetermiWd  by  the  experimenter,  a  value  of  t^  can 
be  obtained  which  is  optimal  in  some  sense.  For  example,  t^  may  be  chosen 
to  minimize  the  AST  when  6  ~ subject  to  constraints  (3.1). 

If  our  purpose  were  to  obtain  a  decision  as*soon  as  possible  without 

A 

regard  to  estimation,  then  other  sequential  tests  based  on  6  could  be 

IB 

constructed  to  meet  this  requirement. 
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